Nagel Scaling, Relaxation and Universality in the Kinetic Ising Model on an 

Alternating Isotopic Chain 
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The dynamic critical exponent and the frequency and wave-vector dependent susceptibility of the 
kinetic Ising model on an alternating isotopic chain with Glauber dynamics are examined. The 
analysis provides to our knowledge the first connection between a microscopic model and the Nagel 
scahng curve originally proposed to describe dielectric susceptibility measurements of several glass- 
forming liquids. While support is given to the hypothesis relating the Nagel scaling to multiple 
relaxation processes, it is also found that the scaling function may exhibit plateau regions and does 
not hold for all temperatures. 



Experimental work on dielectric relaxation in glass- 
forming liquids has in recent years been reported in terms 
of a new (thought to be universal) scaling function |^] 
which is presumed to be related to multifractal scaling. 
While the more usual normalized Debye scaling in terms 
of a single relaxation time is very simple (one chooses 
to scale the frequency with the inverse of the relaxation 
time and the real and imaginary parts are then divided by 
their values at zero and one, respectively), in the so-called 
Nagel plot the abscissa is (1 -|- W) log^Q (to/ujp) jW^ and 
the ordinate is log^Q (x"(w)wp/wAx) jW . Here, x" is the 
imaginary part of the susceptibility, W is the full width 
at half maximum of x" , ^ is the frequency and the one 
corresponding to the peak in x"; and Ax = x(0) — Xoo 
is the static susceptibility. Despite its undeniable phe- 
nomenological success, such scaling is not quite well un- 
derstood on a physical basis. The authors of this proposal 
advance the idea that the presence of more than one re- 
laxation process is not alien to this form and thus suggest 
that multifractality such as the one present in theories of 
chaos may well be behind the new scaling. In order to 
gain some insight into the physical origin of the Nagel 
scaling, it seems appropriate to consider simple but well 
established models in which both universal features are 
unquestionable and more than one relaxation mechanism 
is present. A good candidate may be found among kinetic 
Ising models. 

The scaling hypothesis of Halperin and Hohenberg ||] 
relates the time scale t and the correlation length and 
introduces the dynamic critical exponent z. In the case of 
Ising models, z was for a long time believed to be univer- 



sal, depending on the nature of conserved quantities and 
of those features, for instance dimensionality 1^, which 
determine their static universality class. However, it is 
now well established that for some simple systems this 
exponent is non universal [Q- pO| ]. In particular, in the 
case of one-dimensional Glauber dynamics |l^ the alter- 
nating isotopic chain ||^ presents universal behavior (in 
the sense that it leads to the same value of the dynamic 
critical exponent as the homogeneous chain) whereas the 
alternating-bond chain does not ||l9| (see however 
Ref. ||2^). Due to the fact of this universality of z, the 
isotopic alternating Ising chain with Glauber dynamics 
provides a test model in which to assess the value of the 
connection between multiple relaxation mechanisms and 
the Nagel plot. The model consists of a closed linear 
chain with N sites occupied by two isotopes (character- 
ized by two different spin relaxation times) that are al- 
ternately arranged. The Hamiltonian is the usual Ising 
Hamiltonian given by 



N 



H 



(1) 



where (jj is a stochastic (time-dependent) spin variable 
assuming the values ±1 and J the coupling constant. The 
configuration of the chain is specified by the set of values 
{ci, (72, ■••O'Af} = {c"'^} at time t. This configuration 
evolves in time due to interactions with a heat bath. We 
assume for this chain the usual Glauber dynamics so that 
the transition probabilities are given by 
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(2) 



where 7 = tanh (2 J/fcsT), fcs being the Boltzmann con- 
stant and T the absolute temperature, and is the in- 
verse of the relaxation time of spin i in the absence of 
spin interactions. 

If we now let ai and a2 represent the inverses of the 
free spin relaxation times of chains composed solely of 
spins of species 1 or species 2, respectively, then we can 
set tti = oil — (— l)*a2 , where ai = (ai + a2)/2 and 
7x2 = (ai — Q;2)/2. 

The time dependent probability P(|(t^},<) for a 
given spin configuration satisfies the master equation 



dc^j {t',t' + t) 



dt 



n,j{t')S{t)-ajC,j {t',t' + t) 



+ ^(c,j_i {f, t' + t)+ c,j+i (t', <' + t)). (8) 

We now impose translational invariance and introduce 
Qk , the (spatial) Fourier transform of qj , the t' —f 00 limit 
of the (temporal) Fourier transform of ci{t',t' + t) = 
Cij (t' ,t' + 1) (with I = j — i) denoted by ci{uj), and 
Cfe(w), the spatial Fourier transform of ci{uj), defined 
through 



1j 



(9) 



TV 

+ ^«;,(-a,)P(T,{a^},i), (3) 



where Tijcr^} = {cti, (T2, ...CTi-i, -ctj, (Ti+i, ...o-at}. The 
dynamical properties we are interested in, namely the 
dynamic critical exponent and the susceptibility, re- 
quire the knowledge of some moments of the probability 
P {{cr^^ ,t). Hence, we introduce the following expec- 
tation values and correlation functions defined as: 



(t) = ( 



{<t"} 



(4) 



(5) 



and 



{t', t' + t) = e (t) ((T, {t') a J {f + t)) 



- ^ a^P(K'},tOa,p(K}lK'},0, (6) 

where (t) is the Heaviside step function and the 
sums run over all possible configurations. The sec- 
ond equality of Eq. (|^), which gives the formal defi- 
nition of the time-delayed correlation function, involves 
p (|(T^} I {ct^'I , , the conditional probability of the 
chain having the configuration {o"'^} at time t' ^ t pro- 
vided it had the configuration {cr^'} = {o-^, CTji •■•7 ctat} 
at time t' . Multiplying the master equation by the ap- 
propriate quantities and performing the required summa- 
tions we obtain the set of time evolution equations that 
will be used in our later development. These are given 

by 



dt 



- 2^13-^ 



(7) 



and 



c/(w)= lim — / ci {t' ,t' + t) exp{-iujt)dt (10) 

i'— >oo 27r 



and 



Cfc(w) = (cr-feCTfe)^ = -^^Q(w)exp(-iA:Z). (11) 

In terms of these quantities, Eqs. (^ and (8) may be 
rewritten, respectively, as 



dt 



(12) 



and 



iujci{uj) = rf° - ai {ci{u)) - ^(c;_i(tj) -I- q+i(w))^ 



where 



(13) 



(14) 



_/-«i(l-7C0sfc) -a2(l + 7Cosm 

y —a2 [1 — -f cos k) — ai (1 -f 7Cos/c) J' ^ ' 

and rf° = limt^oo''/(i) is the value of the pair correlation 
function corresponding to the stationary solution of the 
equations of motion in the limit i — > 00. 

The solution to Eq. (|l2|), which yields the magnetiza- 
tion, is straightforward, namely 



(t) 



-Mkt 



*fc (0) 



(16) 



The relaxation process of the wave-vector dependent 
magnetization is determined by the eigenvalues of . 
These are given by 



= -ai±^Ja2 + {al-aflj'^cos'^k. (17) 

The inverses of the (fc-dependent ) relaxation times 
of the ±/c*'* modes are precisely the A^. In the critical 
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region, that is when T and A; ^ 0, A^, —2ai 
while 0. This means that the critical mode is the 

one corresponding to A^. As for the relaxation time, in 
this limit one gets 



Re(-A+)==- 



1 

Tk 



(18) 



2ai 



where we have identified the correlation length ^ as 
^ = exp{2J/kBT)/2. By comparing the former ex- 
pression with the one of the dynamic scaling hypothe- 
sis TjT^ ~ ^~^f{S,k), one finds z — 2, so that, as stated 
above, in the case of one-dimensional Glauber dynamics 
the alternating isotopic chain |^ leads to the same value 
of z as the homogeneous chain. 

Now we turn to the calculation of the other interesting 
response function, namely the frequency and wave- vector 
dependent susceptibility Sk (w), which, by virtue of the 
fluctuation-dissipation theorem [pi|, is defined by 



SkH 



knT 



knT 
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where {<7k<7^k)oo — 1/ (1 — 7C0sfc) cosh(2J//cBr) is the 
static correlation function and {(JkO'^k)^ is the Fourier 
transform of the dynamic one. It should be noted that 
X = ksTSQiui)/ (o'oO'o)oo- After some rather lengthy 
but not too complicated algebraic manipulations start- 
ing with Eq.(|l^) one may arrive at the following result 
P|, namely 



Sk {lo) 



1 



1 - 



ksT (1 — 7 cos k) cosh 

iuj [iuj + ai(l + 7 cos k)) 



{iu! + ai) — ^^'^aia2 (1 -|- cos2A:) — a 



(20) 



from which, using also Eq. ([l7| ) with k = 0, x can be 
expressed in the form 



(1 - j){ai + 02) 



1 -/(«!, 02,7) 1 + /("l, 0:2,7) 



-A^ 



An 



(21) 



Here, the (temperature dependent) function /(ai,Q;2,7) 
is given by 



/(ai,a2,7) 



(ofi — 02)^ — 4aiQ;27^ 



(qi + a2)y/{ai - a-i)^ + 4aia27^ 



(22) 



If we set a\ = a2 in Eq. (21). i.e. we take the uni- 
form chain, then of course the resulting susceptibility has 
the simple Debye form. Although not shown, we have 
checked that this form does not lead to Nagel scaling. 
On the other hand, for the case 7 = 0, we get 



_ 1 

X7=0 - 2 



a2 



ai iuj + 0.2 



(23) 



so that the general structure of the result for the sus- 
ceptibility of the alternating isotopic chain is preserved 
irrespective of the value of 7 (i.e. of the temperature), 
namely a linear combination of two Debye-like terms. 

In Figs. 1 to 3 we present Nagel plots for the cases 
ai = 1 and 02 = 2, ai = 1 and a2 ~ 10, and a\ = \ 
and a2 = 1000, respectively, and different values of 
1/T* = 2J/kBT. In the inserts we include the plots 
x"(u;)/x"(wp) against Lo/ojp which are the natural vari- 
ables of the Debye scaling. We note that as soon as the 
relaxation times become different, except for very low val- 
ues of T* , the agreement with the Nagel scaling improves 
significantly (c/. Fig. 1) until such scaling is virtually 
perfect as depicted in Figs. 2 and 3. On the other hand, 
the almost perfect Debye scaling in Fig. 1 is completely 
lost in Fig. 3. It should be noted that if the two time 
scales are very different, plateau regions eventually ap- 
pear in the Nagel plot as clearly seen in Fig. 3. Whether 
the presence of more than two relaxation times, even if 
not as widely separated as in the case of Fig. 3, would 
also lead to the same type of results is something requir- 
ing future assessment. Also, the precise location of the 
critical value of T* above which the Nagel scaling holds 
as well as the nature of the crossover and of the 'low T*' 
regime are worth investigating. 

It is important to point out that the experiments in 
which the Nagel plots have been more successful concern 
glass-forming systems in which topological constraints 
are assumed to be crucial. However, a clearcut connec- 
tion between such constraints and the different relaxation 
mechanisms has not been established. In this sense, it is 
rewarding that the alternating isotopic Ising chain, which 
is relatively simple with regards to relaxation phenom- 
ena but shows nevertheless universal behavior in terms 
of the dynamic critical exponent, provides perhaps the 
first microscopic model in which this scaling is shown to 
arise. We further want to mention that in order to in- 
clude some of the features present in systems such as 
the above mentioned glass-forming liquids, we have also 
considered in the present context a generalization to an 
. alternating isotopic chain of our quasi one-dimensional 
kinetic Ising- like model of linear polymeric chains 1 22 1 , in 
which the Hamiltonian was chosen as to reduce to the one 
giving the intramolecular energy of the Gibbs-di Marzio 
lattice model Interestingly enough in this model, 

in which the stochastic dynamics implied a rule of tran- 
sition for the configurational changes which was tied to 
the creation or disappearance of flexes and so only some 
states were selected (in the magnetic language this means 
that the domain wall motion is through a biased random 
walk), similar conclusions concerning the Nagel scaling 
readily follow. These will be reported elsewhere. Finally, 
it would be interesting to test whether the same kind of 
scaling is present in other Ising models related to glassy 
systems, such as the spin facilitated kinetic Ising model 
originally introduced by Fredrickson and Andersen p4| ] 
and recently studied in connection with glassy dynamics 
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Figure captions 

Figure 1. Nagel plot for ai = 1 and ai = 2 and for 
T* = 1,2,5, 10 and 100. There is reasonable agreement 
with the scaling form for this choice except for low T* . 
The insert contains the plot of x"(cL')/x"(a;p) vs. w/wp 
in order to test the Debye-like behavior. 

Figure 2. The same as Fig. 1 but with the choice 
ai = 1 , = 10 and T* = 5, 10, 50 and 100. The 
improvement in the agreement with the Nagel scaling is 
rather noticeable, while the opposite trend is observed 
with respect to the Debye scaling. 

Figure 3. The same as Figs. 1 and 2 but for ai = 1 , 
a2 = 1000 and T* = 5, 10, 50 and 100. A plateau region 
in the Nagel plot is clearly present in this case. Here, the 
behavior is definitely non-Debye. 
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